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The study of diboson production in proton-proton collisions provides an important test of the 
standard model (SM). Many extensions of the SM predict new scalar, vector, or spin-2 particles 
that decay into a pair of W or Z bosons. In addition, these final states are sensitive to the self- 
interactions among the gauge bosons via trilinear gauge couplings (TGCs). These couplings 
are the direct consequence of the non-Abelian SIT (2) x LT(1) gauge symmetry of the SM and 
are a necessary ingredient to construct renormalizable theories. The values of these couplings 
are fully determined in the SM by the gauge structure of the Lagrangian. Therefore, any devi- 
ation of the observed coupling strength from the SM prediction would indicate the presence of 
new physics. This deviation would be manifested as a change in the production cross section, 
especially for energetic heavy gauge bosons. 

In the SM, ZZ production proceeds via the t- and u-channel qq scattering diagrams, and via 
gluon-gluon fusion. The presence of anomalous neutral trilinear couplings (ATGCs) would 
lead to a sizable enhancement of ZZ final states via s-channel qq scattering. A model featuring 
such couplings can be constructed by means of an effective Lagrangian 111. In this parametriza- 
tion, two ZZZ couplings and two ZZ7 couplings are allowed by electromagnetic gauge invari- 
ance and Lorentz invariance for on-shell Z bosons. The couplings are parametrized by two 
CP-violating (ff'^) and two CP-conserving iff' 7 ) complex parameters, which are zero in the 
SM. 

Measurements of the ZZ cross section were previously performed at the Tevatron El [3j and 
the Large Hadron Collider (LHC) |@H5|. A first measurement of the ZZ cross section at center- 
of-mass energy y/s = 7 TeV by the Compact Muon Solenoid (CMS) Collaboration in the decay 
mode ZZ — > 2£2£, where £ is either e or ji, is presented in Ref. |U. The measured cross section 
£7(pp -»■ ZZ)£(ZZ -> 2£2£) = 28.1 j^g (stat.) ± 1.2 (syst.) ± 1.3 (lumi.) fb agrees well with the 
SM prediction of 27.9 ± 1.9 fb. In this Letter, we present an extended measurement of the ZZ 
production cross section based on the decay mode 2£2£' , where £' is e, ji, or t. If a T is present in 
the final state, one Z is required to decay into e + e~ or , and the second Z into t + t~ in four 
possible final states: ThTh, T e Th, t^t^, and t 6 t^, where Th represents a x decaying hadronically 
while T e and T f , indicate taus decaying into an electron and a muon, respectively. The presence 
of four leptons in the final state provides a clean signature with only a small contribution from 
background processes. The background sources include reducible contributions from Zbb and 
tt processes, where the final states contain two isolated leptons and two b jets with secondary 
leptons, and from Z+jets and ZW+jets processes where the jets are misidentified as leptons. 

In this Letter, we also present a search for the neutral ZZZ and ZZ7 ATGCs. Previous studies 
on neutral ZZZ and ZZ7 ATGCs were performed at LEP2 I6 HT0I , the Tevatron l|TT| , and the 
LHC 0. The most restrictive limits were set in Ref. 0, -0.07 < < 0.07 and -0.08 < 
/7g < 0.08, with a data set corresponding to an integrated luminosity of 1 fb -1 of pp collisions 
atVs = 7 TeV. 

The measurements presented here are based on data collected in 2010 and 2011 with the CMS 
experiment at the LHC at \Js = 7 TeV, corresponding to an integrated luminosity of 5.0 ± 
0.1 fb -1 . A set of Monte Carlo (MC) event samples is used to simulate signal and background 
events. The ZZ production via qq is generated at next-to-leading order (NLO) with POWHEG 
2.0 1(12^1411 , while other diboson processes (WW, WZ, Z7) are generated with pythia 6.424 
and MadGraph ||T5 | . The gg — > ZZ contribution is estimated using events generated with 
the gg2ZZ code JTSJ. The Z+jets samples, namely Zbb , Zee, and Z+light jets, are generated 
with MADGRAPH. The tt events are generated at NLO with POWHEG. For leading-order gener- 
ators, the default set of parton distribution functions (PDFs) used to produce these samples 
is CTEQ6L JT7|, while CT10 Ell is used for NLO generators. Finally, for the modeling of 



2 



ATGCs, the SHERPA generator version 1.2.2 is used [19 j. The r-lepton decays are generated 
with TAUOLA fl20|l . All events are processed through a detailed simulation of the CMS detector 
based on Geant4 Ij2l) and reconstructed with the same algorithms as used for data. 

A detailed description of the CMS detector can be found elsewhere |22| . The central feature of 
the CMS apparatus is a superconducting solenoid of 6 m internal diameter, providing a mag- 
netic field of 3.8 T. Within the field volume are the silicon pixel and strip tracker, the crystal 
electromagnetic calorimeter (ECAL), and the brass /scintillator hadron calorimeter (HCAL). 
Muons are measured in gas-ionization detectors embedded in the steel flux return yoke of the 
magnet. CMS uses a right-handed coordinate system, with the origin at the nominal interaction 
point, the x axis pointing to the center of the LHC ring, the y axis pointing up (perpendicular 
to the plane of the LHC ring), and the z axis along the counterclockwise-beam direction. The 
polar angle 6 is measured from the positive z axis and the azimuthal angle cp is measured in 
the x-y plane. Variables used in this analysis include the pseudorapidity rj = — In [tan (6/2)] 

and the transverse momentum pj = \jp\ + Py- The ECAL is designed to have both excel- 
lent energy resolution and high granularity, properties that are crucial for reconstructing elec- 
trons and photons produced in r-lepton decays. The ECAL is constructed with projective lead 
tungstate crystals that provide coverage in pseudorapidity \tj\ < 1.48 in a barrel region and 
1.48 < | ?/ 1 < 3.00 in two endcap regions (EE). A preshower detector consisting of two planes 
of silicon sensors interleaved with a total of 3Xo of lead is located in front of the EE. The en- 
ergy resolution is 3% or better for the range of electron energies relevant for this analysis. The 
tracker measures charged particle tracks within the range \t]\ < 2.5. It consists of 1440 silicon 
pixel and 15 148 silicon strip detector modules, and provides an impact parameter resolution of 
~15 pm and a transverse momentum resolution of about 1.5% for 100 GeV particles. The recon- 
structed tracks are used to measure the location of interaction vertices. The spatial resolution 
of the reconstruction is ~25 pm for primary vertices with more than 30 associated tracks 11231 . 
The barrel region of the muon system is instrumented with drift tubes, and the endcap regions 
with cathode strip chambers. In both regions, resistive-plate chambers provide additional co- 
ordinate and timing information. Muons are reconstructed in the range \r]\ < 2.4, with a typical 
px resolution of ~1% for pj = 40 GeV. 

At the trigger level, the selected events are required to have either at least two electrons, one 
with pj > 17 GeV and the other with pj > 8 GeV, or at least two muons, one with pj > 13 GeV 
(px > 17 GeV for high instantaneous luminosity data-taking periods) and the other with pj > 
8 GeV. 

Electrons are reconstructed within \t] e \ < 2.5 and with p| > 7 GeV by combining information 
from the ECAL and tracker Il24l |25| . Electron identification requirements rely on the electro- 
magnetic shower shape and other observables based on tracker and calorimeter information. 
The selection criteria depend on p| and \f] e \, and on a categorization according to observables 
that are sensitive to the amount of bremsstrahlung emitted along the trajectory in the tracker. 
Muons are reconstructed |26| within \rjl'\ < 2.4 and p T > 5 GeV with information from both 
the tracker and the muon spectrometer. The track must have more than 10 hits in the silicon 
tracker |22H to ensure a precise measurement of the momentum. The efficiencies are measured 
in data, using a tag-and-probe technique |27] based on an inclusive sample of Z — > £ + £~ events. 
The measurements are performed in several ranges of p T and \rj |. The product of reconstruc- 
tion and identification efficiencies for electrons in the ECAL barrel (endcaps) varies from about 
68% (62%) for the p\ range 7-10 GeV, to 82% (74%) at 10 < p\ < 20 GeV, and reaches up to 90% 
(89%) at p\ > 20 GeV. The muons are reconstructed and identified with efficiencies above 98%. 

Since the ZZ final state is expected to have only a small contribution from background pro- 
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cesses, the algorithms are tuned to maximize the lepton-reconstruction efficiency, resulting in 
an increased lepton-misidentification rate. A particle-flow (PF) technique l|28| is used for Th 
reconstruction. In the PF approach, information from all subdetectors is combined to recon- 
struct and identify particles produced in the collision. The particles are classified into mutually 
exclusive categories: charged hadrons, photons, neutral hadrons, muons, and electrons. These 
particles are used to reconstruct the Th candidates with the "hadron plus strip" (HPS) algo- 
rithm [29 1, which is designed to optimize the performance of Th identification and reconstruc- 
tion by considering specific Th decay modes. The neutrinos produced in all x decays escape 
detection and are ignored in the Th reconstruction. The algorithm provides high Th identifica- 
tion efficiency, approximately 50% for the range of Th energies relevant for this analysis, while 
keeping the misidentification rate for jets at the level of 1%. 

Events are required to have at least one Z — > l + £~ candidate, denoted by Zi. The invariant 
mass of the reconstructed Z\ is required to be 60 < tnu < 120 GeV. The two leptons must have 
opposite charges, one with pj > 20 GeV and the other with pj > 10 GeV, and with \ t]\ < 2.5 for 
the electrons and \rj\ < 2.4 for the muons. If more than one candidate is found, the one with 
the mass closest to the Z mass is considered as Z\. 

Lepton isolation requirements depend on the ZZ decay mode. For the final states with only 
electrons and muons, the isolation criteria are based on a combination of the tracker, ECAL, 
and HCAL information. The standard combined relative isolation is defined as 



rstd 
1 rel 



E PT, track + m ax(E 4, ECAL + E E X HCAL " 71 " AR max * P) 0) J /pi 



with the sums running over the charged tracks and the energy deposits in the ECAL and 
HCAL within a cone around the lepton direction defined by AR < AR max = 0.3, where 
AR = y/ (At]) 2 + (A</>) 2 , and Ej stands for the transverse energy. The isolation is made largely 
insensitive to pileup of pp collisions by correcting for the average energy density, p, measured 
in cones at randomly chosen locations to sample the average underlying energy in the event. 
The isolation variable I^f is required to be less than 0.275 for each lepton. The significance 
of the impact parameter of each lepton relative to the event vertex (Ssd) is required to satisfy 
I S3D I < 4. The primary vertex is chosen as the vertex with the highest sum of p\ of its con- 
stituent tracks. 

In the 1l2x final states, instead of standard isolation, the leptons from the Zi are required to 
have a combined PF relative isolation 1^ < 0.25. The 1^ is defined similarly to I^f, however 
in this case the sums run over charged hadrons, photons, and neutral hadrons, all measured in 
the isolation cone of AR < 0.4 around the lepton direction. 

The selection requirements for the second Z, denoted by Z2, also depend on the final state. In 
the final states with electrons and muons only, the isolation requirements are the same as for 
the leptons from Z\, but pj > 7 GeV and pj > 5 GeV are required for electrons and muons, 
respectively. If the final state is t 6 t^, the lepton pj values are required to exceed 10 GeV. The 
remaining criteria are identical to those for Z\. Since hadronically decaying x leptons have 
much larger misidentification rates than the other leptons, the isolation requirement based on 
^rei f° r * ne electrons and muons in the final states T e Th and T^Tj, is changed to 0.15 and 0.1, 
respectively. A study of inclusive Z — > t + t~ production |l30l| demonstrated that modifying 
the electron and muon isolation requirements is a more effective way to reduce background 
in such final states than requiring tighter isolation on Th- The x leptons are required to have 
pj > 20 GeV and \rj\ < 2.3, and to satisfy the requirements of a loose HPS working point. If 
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the Z2 decays to , both Th are required to satisfy the requirements of a medium working 
point of the HPS algorithm. The loose (medium) working point requires the scalar sum of 
the px of the charged hadrons and Ej of the neutral hadrons within the isolation cone to be 
less than 2GeV (lGeV). The loose (medium) working point corresponds to a probability of 
approximately 1% (0.5%) for jets to be misidentified as t^. Using the medium instead of loose 
working point leads to a decrease in the Th reconstruction efficiency from ^50% to «40%. 

The invariant mass of the reconstructed Z2 is required to satisfy 60 < mn < 120 GeV, when 
Z2 decays into e + e~ or p + p ~ . In the 2£2t final states, the visible invariant mass of the recon- 
structed Z2 — > t + t~ is required to satisfy 30 < m" T s < 80 GeV. The upper bound reduces 
contributions from Z2 — > £ + £~ , where an electron or a muon is not well reconstructed and is 
misidentified as a Th. For the Z2 — > t 6 t^ final state, the upper bound on m TT is increased to 
90 GeV, as this state is not produced in Z2 — > £ + £~ decays. In the final states involving Th, 
leptons from the same Z are required to be separated by AR > 0.4 for the Zi, and by AR > 0.5 
for the Z2. 

The major contributions to the background are due to Z production in association with jets, 
WZ production in association with jets, and tf . In all of these cases, a jet or nonisolated lepton 
is misidentified as an isolated electron, muon, or Th- The relative contribution of each source of 
background depends on the final state. 

The background estimate is performed in two steps. Firstly, the rate for loosely isolated objects 
to be misidentified as isolated ones is measured in a control region that does not contain any 
signal contribution. The misidentification rate is estimated with events in which the Z\ passes 
all selection requirements, and which contain an additional probe electron, muon, or Th- No 
isolation requirement is applied to the probe. The misidentification rate is defined as the ratio 
of the number of probe candidates that pass the isolation requirements to the initial number of 
probe candidates, and is measured as a function of px and rj for each lepton flavor. 

The second step is to estimate the number of background events in the signal region. The mea- 
sured misidentification rate is applied to events that pass all selection requirements, including 
the opposite-charge requirement for the Z2, but requiring the candidates to not be isolated. 

Theoretical uncertainties on the ZZ — > 2£2£' acceptance are evaluated using MCFM 6.2 pl|, 
varying QCD scales up and down by a factor of two with respect to the default factorization 
(ji?) and renormalization (^r) scales pp = Pr = mz, where mz is the mass of the Z boson. The 
variations in the acceptance are 0.1% (qq — > ZZ) and 0.4% (gg — > ZZ) and can be neglected. 
The uncertainties related to the PDFs are evaluated following the PDF4LHC prescription II32I . 
Using the CT10 HI, MSTW08 |331, and NNPDF |]3U sets, the uncertainties are estimated to be 
4% for qq — > ZZ and and 5% for gg — > ZZ processes. 

The uncertainties on Z+jets, WZ+jets, and tt backgrounds reflect the uncertainties on the mea- 
sured values of the misidentication rates and the limited quantity of data in the control regions 
in the data and amount to 30-50% depending on the decay channel. The uncertainty on the 
integrated luminosity is 2.2% [35J. Systematic uncertainties on trigger efficiency (1.5%), lepton 
identification efficiency, and lepton isolation are evaluated from data. The uncertainties asso- 
ciated with lepton identification and isolation are 1-2% for muons and electrons, and 6-7% for 
Th. Uncertainties on energy scales, 3% for Th and 1-2.5% for electrons, contribute to variations 
in the shape of the mass spectrum. 

Table [l] presents the number of observed events in the signal region in each channel, as well as 
the expected number of signal events and the estimated number of background events. We ob- 
serve a total of 54 candidate events in the 4e, 4p, and 2e2p channels, compared to the SM expec- 
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Table 1: The expected yield of ZZ events obtained from simulation and the estimated yield 
of background events obtained from data, as described in the text, are shown for each decay 
channel and are summed in the total expected yield ("Total expected"). They are compared to 
the number of events observed in the signal region. The first uncertainty is statistical while the 
second one is systematic. 



Decay 


Expected ZZ 


Background 


Total 


Observed 


channel 






expected 




eeee 


10.50 ± 0.04 ± 0.95 


0.25 ±0.14 ±0.07 


10.75 ± 0.14 ± 0.95 


9 


mm 


15.91 ± 0.05 ± 1.43 


0.52 ± 0.26 ± 0.25 


16.43 ±0.26 ±1.45 


14 


eejiji 


26.74 ± 0.10 ± 2.41 


0.58 ±0.18 ±0.23 


27.32 ± 0.21 ± 2.41 


31 


eeT h T h 


0.75 ± 0.01 ± 0.07 


0.76 ±0.16 ±0.05 


1.51 ±0.16 ±0.09 


1 


m r h^h 


0.82 ± 0.02 ± 0.07 


0.75 ±0.16 ±0.08 


1.57 ±0.16 ±0.11 





eeT e T h 


1.17 ±0.02 ±0.11 


0.96 ± 0.34 ± 0.12 


2.29 ±0.34 ±0.16 


3 




1.15 ±0.02 ±0.10 


0.35 ±0.34 ±0.11 


1.60 ±0.34 ±0.15 


3 


eer f( T h 


0.94 ± 0.02 ± 0.08 


0.22 ±0.14 ±0.04 


1.17 ±0.14 ±0.06 





m T \^ 


1.08 ±0.02 ±0.10 


0.55 ±0.24 ±0.11 


1.64 ±0.24 ±0.15 


2 


eer e T fl 


0.54 ± 0.01 ± 0.05 


0.64 ± 0.44 ± 0.16 


1.22 ±0.44 ±0.17 





m T ^Y- 


0.60 ± 0.01 ± 0.05 


0.14 ±0.30 ±0.10 


0.74 ±0.30 ±0.11 


2 



tation of 54.5 ± 0.3 (stat.) ± 4.8 (syst.) events, which includes 1.4 from background processes. In 
the 2£2t channels, 11 candidate events are observed, compared to 11.7 ± 0.8 (stat.) ± 1.0 (syst.) 
events expected, including 4.4 from background processes. The reconstructed four-lepton in- 
variant mass distributions are compared to the SM expectations in Figs, [l] (a) and (b) for the 
sum of the 4e, 4^/, and 2e2}i channels, and the sum of all the 2£2t channels. The shapes of the 
signal and background are taken from the MC simulation, with each component normalized 
to the corresponding estimated value from Table [lj The reconstructed masses in 2£2t states 
{m™2 T ) are shifted downwards with respect to the generated Z masses by about 30% due to the 
undetected neutrinos in r decays. Figures [l] (c) and (d) demonstrate the relationship between 
the reconstructed Zi and Z2 masses. 

To include all the final states in the calculation of the cross section, a simultaneous fit to the 
numbers of observed events in all the decay channels is performed. The fit is constrained by 
the requirement that all the measurements come from the same initial state via different de- 
cay modes. It allows for combining many decay modes with either very few or no events 
observed. The joint likelihood is a combination of the likelihoods for the individual channels, 
which include the signal and background hypotheses. Each r-lepton decay mode is treated as 
a separate channel because they are mutually exclusive owing to the methodology adopted for 
event reconstruction and subsequent event selection. The statistical and systematic uncertain- 
ties are introduced in the form of nuisance parameters via log-normal distributions around the 
estimated central values. 

The resulting cross section is measured to be 

<7-(pp -> ZZ) = 6.24 ±g-fg (stat.) t° H (syst.) ± 0.14 (lumi.) pb. 

This result is to be compared to the theoretical value of 6.3 ± 0.4 pb calculated with MCFM at 
NLO for qq -> ZZ and LO for gg -> ZZ with the MSTW08 PDFs and for both Z bosons in the 
mass range 60 < Mz < 120 GeV. This is the most precise published pp — » ZZ cross section 
measurement to date, which for the first time extends the pp — > ZZ — > 2£2£ measurement to 
include final states with hadronically decaying r leptons. 
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Figure 1: Distributions of the four-lepton reconstructed mass for (a) the sum of the 4e,4^, and 
2e2}i channels and (b) the sum of the 212t channels. Points represent the data, and the shaded 
histograms represent the expected ZZ signal and the reducible background. The shapes of the 
signal and background are taken from the MC simulation, with each component normalized to 
the corresponding estimated value from Table [Tj The distributions (c) and (d) demonstrate the 
relationship between the reconstructed Zi and Z2 masses. 
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Figure 2: Distribution of the four-lepton reconstructed mass for the sum of the 4e,4^, and the 
2e2fi channels. Points represent the data, and the shaded histograms represent the expected ZZ 
signal and the reducible background. The dashed and dotted histograms represent the results 
of the SHERPA simulation for the SM (/f = 0) and in the presence of an ATGC (/f = 0.015), 
while all the other anomalous couplings are set to zero. 



CMS 



0.04 
0.02 



-0.02 



-0.04 



= 7 TeV, L = 5.0 fb" 1 

Observed pp^zz->2I2I _ 
Expected ± 1a 95% CL 
Expected ± 2a 

fc^ (a) - 




No form factor assumed 
ATGC values outside contour excluded 



-0.02 -0.01 



0.01 0.02 



CMS 



0.04 
0.02 



-0.02 



-0.04 



\[s = 7 TeV, L = 5.0 fb" 1 

Observed ' PP^zz-»2i2i 
Expected ± 1a 95 . /o CL 
Expected ± 2a 




No form factor assumed 
ATGC values outside contour excluded 

J i l i i 



-0.02 -0.01 







0.01 0.02 



Figure 3: Expected and observed two-dimensional exclusion limits at 95% CL on the anomalous 
neutral trilinear ZZZ (/f 5 ) and ZZ7 (fj 5 ) couplings. The green and yellow bands represent the 
one and two standard-deviation variations from the expected limit. In calculating the limits, 
the anomalous couplings that are not shown in the figure are set to zero. 
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The limits on ATGCs are calculated with the modified frequentist construction CL S |36ti38| 
based on the shape of the four-lepton invariant mass distributions, including the 4e, 4^, and 
2e2ji channels in the likelihood combination. Figure [2] presents the distribution of the four- 
lepton reconstructed mass for the sum of the 4e, Aji, and 2e2ji channels. The dashed and dotted 
histograms represent the results of the SHERPA simulation for the SM (ff = 0) and in the 
presence of an ATGC (ff = 0.015), while all the other anomalous couplings are set to zero. The 
presence of ATGCs would be manifested in an increased yield of events at high four-lepton 
masses. The invariant mass distributions are interpolated from SHERPA simulation for different 
values of the anomalous couplings. For each distribution only one or two couplings are varied, 
while all others are set to zero. The fit is performed to find the maximum likelihood value and 
limits are calculated. To avoid unitarity violation at energies above the scale A of new physics, 
the ATGCs are often modified with a form-factor parametrization of the type 1/ (1 + s/A 2 ) 2 , 
where \/s ~ 1112121 is the effective center-of-mass energy of the collision. However, no unitarity 
violations occur in the sensitive region m^ni ^ 1.5 TeV for bare anomalous couplings of order 
0.05 or smaller II39II , so we calculate the limits without form-factor scaling. This choice has the 
advantage of avoiding any bias from energy-dependence assumptions and is exact in the limit 
in which the scale of new physics is much larger than y/§. 

Figure [3] presents the expected and observed two-dimensional exclusion limits at 95% confi- 
dence level (CL) on the anomalous neutral trilinear ZZZ and ZZ7 couplings. The green and 
yellow bands represent the one and two standard-deviation variations from the expected limit. 
The present limits are dominated by statistical uncertainties. Systematic uncertainties arising 
from the uncertainty on the theoretical cross section, PDFs, detector efficiencies, and luminosity 
are introduced in the form of nuisance parameters with log-normal probability density func- 
tions. One-dimensional 95% CL limits for the ff ,ry and /i?' 7 anomalous coupling parameters 
are measured to be 

- 0.011 < ff < 0.012, -0.012 < ff < 0.012, -0.013 < fj < 0.015, -0.014 < / 5 7 < 0.014. 

In the one-dimensional fits, all of the ATGC parameters except the one under study are kept 
fixed to zero. These values extend previous results on vector boson self-interactions and are 
currently the most stringent limits established for ZZZ and ZZ7 couplings. 

In summary, we have presented an updated measurement of the ZZ production cross section 
in proton-proton collisions at 7 TeV in the ZZ — > 2£2£' decay mode, with £ = e, p and £' = 
e,fi,r. The data sample corresponds to an integrated luminosity of 5.0 fb -1 . The measured 
cross section cr(pp ->• ZZ) = 6.24 (stat.) t° H (syst.) ± 0.14 (lumi.) pb is consistent with the 
SM prediction and is the most precise published pp — > ZZ cross section measurement to date. 
For the first time the pp — > ZZ — > 2121 measurements are extended to include final states 
with hadronically decaying x leptons. Limits on vector-boson self-interactions are established, 
significantly restricting anomalous ZZZ and ZZ7 trilinear gauge couplings. 
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